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1. INTRODUCTION 

Device to device (D2D) communication has recently been introduced in 3rd-generation partnership 
project (3GPP) long term evolution (LTE) Release 12 and 13. D2D, together with the core technologies, is 
proposed for forthcoming 5th generation (5G) cellular standard. In wireless cellular networks, the device-to- 
device (D2D) communications allow devices in close proximity to transmit and receive signals directly without 
connecting with base stations (BSs). Therefore, to reduce the high load from core networks, D2D has been 
recognized as a reasonable technique [1]-[3]. Together operation of the cellular network, the D2D devices 
are able to share the same spectrum resources with cellular users (CUs), and the spectral efficiency can be 
enhanced. In principle, D2D and cellular communication links are categorized into two kinds, including the 
overlay and underlay spectrum sharing techniques [4]. D2D users can only utilize the idle spectrum that is not 
currently employed by the CUs in the overlay D2D mode. The underlay spectrum sharing mode permits D2D 
pair users to occupy the resource of the cellular connections. The co-channel interference between the D2D and 
cellular links occurs, and such limitations need to be addressed. Therefore, the resource allocation problems 
are studied in D2D communications underlaying the cellular networks [5]-[7]. In [5], to maximize the sum rate 
of D2D links, the authors presented a power allocation scheme, and guarantee the performance of the cellular 
link at the same time. In order to maximize the rate of the single D2D link in [6], power allocation for D2D 
and cellular links is studied. The multiple cellular and D2D links are investigated in [7]. 

Considering the advantages of D2D communication, There are some studies on the usage of D2D com- 
munications in [8]-[12]. The cellular networks require high maintenance costs and additional infrastructure, 
and more difficulties with a large amount of equipment. In order to address these challenges, D2D communica- 
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tions have been proposed as an essential technique to implement with non-orthogonal multiple access (NOMA) 
networks, the heterogeneous network recently [13]-[19]. 

To study a large number of nodes at various locations in D2D, stochastic geometry is a sophisticated 
tool from mathematics to provide insights into spatial averages taken over several realizations. Many metrics 
can be evaluated, such as outage probability, interference, data rate, and signal-to-interference-plus-noise ratio 
(SINR) [20],[21]. Stochastic geometry can be employed to design many applications related to D2D [22],[23]. 
Moreover, the published works [24]-[28] also study the performance of underlay D2D networks. The outage 
probability of random D2D communication in millimeter-wave cellular networks is studied in [24]. Moreover, 
the authors in [25] investigate the performance over nakagami-m fading channels. The paper [26] also attempts 
to investigate the secure property of the D2D network. An interesting model where multiple D2D pairs are 
deployed is proposed in [27]. An extension to full-duplex communication for the D2D network is introduced 
by [28]. In the above works, there is no concern over the spatial separation between the cellular network and 
the D2D network. In addition, statistical properties of the distance between users in cellular networks and D2D 
networks under the nodal distancing constraint have not been reported in the literature. 

Motivated by these analysis, this paper examines outage performance of traditional users located in 
outside of D2D serving area in cellular network. 


2. SYSTEM MODEL 

We consider a downlink wireless network in which there is a base station (BS), a cellular user (CUE) 
sharing resources with a D2D transmitter (D2D-Tx), as plotted in Figure 1. The CUE is uniformly distributed 
inside the cellular coverage area with radius Re around the BS. In addition, the CUE is located outside the D2D 
coverage area with radius Ra < Re around the D2D-Tx. We facilitate each node with a single antenna. 

Further, the channel coefficients of the links from the BS and from the D2D-Tx to the CUE are de- 
noted by h and g, respectively. We assume all channels to be experiencing independent Rayleigh fading with 
‘[|h|?] = E[|g|?] = 1. To capture the large-scale propagation, we adopt the power-law path-loss model where 
R and L are the corresponding distances to the CUE, and 7 denotes the pathloss exponent. Moreover, we 
assume a block-fading model in which all channels are invariant during the transmission of a block and vary 
independently among different blocks. 
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Figure 1. The proposed underlay D2D communication network. 


In the downlink, the CUE receives the information signal from the BS while being interfered by the 
D2D-Tx. Subsequently, the received signal at the CUE is given by 
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in which P, and P, are the transmit powers from the BS and the D2D-Tx, respectively, x and s are the corre- 


























sponding signals. Herein, E[|x|?] = E[|s|?] = 1, and n denotes the additive white Gaussian noise (AWGB) at 
the CUE with zero mean and variance a. Hence, the achievable instantaneous rate at the given CUE is 
p, RÊ 
R = log, 1+ , (2) 
P, Ir + o? 


3. DISTANCE DISTRIBUTION 

This section studies the stochastic properties of the link distances. Particularly, the CDF and PDF of 
the distances R and L are derived in closed-form. It should be mentioned that both L and R are dependent on 
the location of the CUE. 


3.1. Preliminaries 

In this subsection, we specify the distribution, i.e., CDF and PDF, of the distance from a fixed point A 
to the random point C'. Particularly, the point C uniformly distributed inside the circular area with the radius R 
around a fixed point B, i.e., C € C(B, R). The distance probability distribution function is considered in two 
different cases as depicted in [29, Eq.(34)], and can be formulated as 

IC(A, AC) N C(B, R)| 
Fac(AC) = , AB— R< AC < AB +R, (3) 
IC(B, R)| 

in which | - | denotes the Lebesgue measure. Denoting Aror = C(A, AC) N C(B, R), for the case AB > R, 
|Aror| is calculated as by [29] 











2 2_ p2 Der 2_ p2 
JAror| =5(AC? + R?) — AC? sin (22 ee Z) + Resin? (= an =) 


2AC.AB 2R.AB 
— 5 VAC? — (AB — R)?)((AB + R)? — AC?) £ S(AB, R; AC), (4) 





with AB — R < AC < AB + R. For the case AB < R, the area of Aror is obtained as [29] 


(5) 


vera mn AC? ,0<AC<R-AB 
TOT! |) S(AB,R; AC) , R-AB< AC<R+AB 


3.2. The information links 

Initially, let Ac be the 2D-Euclidean space in which the CUE uniformly distributed and Ap be the 
coverage area around the D2D-Tx. To reduce the cross-mode interference, the interference from the D2D-Tx 
to the CUE and vice versa, it is assumed that Ac and Ap are disjoint meaning that the CUE is located outside 
the D2D proximity. Hence, the CDF of R is given as 


_ |C(BS,r)\Ap| 


PRO = Tel 





,0<r< Re. (6) 


Denoting S = C(BS, r)\ Ap, the above expression can further be derived via the four following cases. 
a) Case Dp + Ra < Re ^ Di > Ra: 

As shown in Figure 2, |Ac| = 7R? — tR? and the calculation for |S| is divided into three following 
subcases: 1) if 0 < r < Di — Ra, it is clearly that |S| = mr?, 2) if Dr + Ra < r < Re, |S| = nr? — wR? and 
3) if Dp — Ra < r < Di + Ra, thus 


|S| = |C(BS,r)| — |C(BS,r) NAp| = mr? — S(D;, Ra;r). (7) 


Now, the corresponding PDF is obtained by substituting (7) into (6) and then taking the derivative of 
the resulted F’p(r), thus 


2rr „0 <r < D- Ra 
ae Deane (Ca) De Ra<r< D,+ Ra (8) 


2Dir 
, Di + Ra<r < Re 
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AB? + AC? — R? 
) (9) 


2AB.R 


b) Case Di + Ra < Re A Di < Ra: 

As shown in Figure 3, |Ac| = mR? — mR? and |S| is derived as following subcases: 1) if 0 < r < 
Ra ~ Dr, thus |S| = 0, 2) if Di + Ra < r < Re, thus |S| = nr? — wR}, and 3) if Dy — Ra < r < Di + Ra, 
|S| = ar? — S(D,, Ra;r). 

As a result, the corresponding PDF is given by 


0 ,0<r<Ra-D; 
1 ide 
fr(r) = aR nE mr +2rsin7! (Poa) , Di — Ra < r < Di + Ra (10) 
° d | orr Dirt Re ere Re 


c) Case Di + Ra > Re > Di — Ra: 
The cellular coverage is given as |Ac| = mR? — |A| where A = C (BS, R.) O Ap as shown in Figure 
4. It is known that [30] 


D? + R} — R? D? + R? — R? 
L p2 anal [ Yt d— fe Paget | Yt c d 
|A| =R4 cos ( 2D, Ra ) + Ré cos ( ZDR, ) 





1 
= ND + Ra + Re)(Di + Ra — Re) (Di — Ra + Re)(Di + Ra + Re). (11) 


In this case, |S| can be derived via ne two following subcases: 1) if0 < r < Ra— Drs, thus |S| = mr’, 
and 2) if Di — Ra < r < Di + Ra, |S| = mr? — S(D;, Ra;r). 
Hence, the PDF of the distance of the information link is then obtained as 


1 2nr „0 <r < Dı- Ra j 
fr(r) = TR2— |A| \ rr + 2rsin7! (P) , Di— Ra<r < Re oal 
d) Case Di — Ra > Re: 

In this case, C(BS, Re) MN Ap = 0, thus the CDF and the PDF of the distance R are immediately given 








as 
r2 
Fr(r) = Re’ 0<r< Re, (13) 
2r 
frír) = RZ’ 0<r< Re. (14) 
7 i / 
/ Ti \ Ra 
< R. 4 È | 
a a 
N oe | < 
fa 7, A realization of Figure 3. A realization of Figure 4. A realization of 
Di + Ra < Re ^ Di > Ra Di + Ra < Re ^ Di < Ra Di + Ra = Re = Di — Ra 
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3.3. The interference links 
In this subsection, the CDF and PDF of the link distance from the interfering D2D-Tx to the CUE are 
studied. Similar to the previous parts, the CDF of the distance L is given as 





_ |AgNC(D2D-Tx, r)| 


Fr(r) |Ac| 


,0<r<D4+ Ra- (15) 


Note that we reuse the notation S ê Ac  C(D2D-Tx, r) for the following analysis. Subsequently, 
the above expression can further be derived via the three following cases. 
a) Case Di + Ra < Re: 

Similarly, |Ac| = mR? — wR? and |S] is derived by the following subcases. 1) if r < Ra, it is 
obvious that |S| = 0, 2) if Ra < r < Re — Di, |S| = mr? — wR?, and 3) if Re — Di < r < Re + Dz, 
|S| = S(Dz, Re; r) — mR}. Subsequently, the PDF of the interference distance is then given as 


0 ,0<r< Ra 
Qrr , Ra<r < Re- Di (16) 


Dine a2 
mr — 2r sin! (Pipe te ) , Re— Di <r < Re+ Di 


1 
=- wR? — rR? 


f(r) 


b) Case D; + Ra > Re > D; = Ra: 
Similar to the previous subsection, |Ac| = mR? — |A| and |S| = S(Di, Ra;r) — |A| if Ra < r < 
Ra + D; whereas |S| = 0 whenever r < Ra. Hence, 





1 0 ,0<r< Ra i 
fir) = TR2— |A| \ ar —2rsin7! (i) , Ra<r < Re+ Di (L) 
c) Case Di — Ra > Re: 
In this case, |Ac| = 7R2 and |S| is obtained as [...] 
S(Dz, Re; , Di — Re < r< Dit Re 
jaja d iR] A PAE (18) 
0 5 0 LENE Di iai Re 
Substituting the above result into (15) and with the help of (9), one can the obtain 
1 nr — 2r sin”! (=) ,D,-Re<r< D+ Re 
fi(r) = z t (19) 
TRe (0 „0 <r < Di- R 


4. OUTAGE PERFORMANCE 

In this section, we investigate the outage probability of the proposed system. Let Rin be the target data 
rate of the CUE, thus the probability for the instantaneous rate at this CUE falls below Rin is mathematically 
expressed as 


Pout =P1r{R < Rin} (20) 


nA x 2 
=1— Pr f in? a hk (4 +1}, (21) 
Yo Ln” 





in which yn £ 2° — 1, % £ P,/o? and 7 £ P,/o?. Now that |h| and |g| are both independent and identical 


Rayleigh distributed random variables with unit variance, thus the CDF of |h|? and the PDF of |g|? are given 
by 


Finj2 (2) =l-e*,2«>0, 
fige (2) =e, > 0, 
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respectively. Hence (21) can further be expressed as 


a Rn (ll 
Pout =] — “R,|g|2,L lew { Ap ( La + 1 (22) 


c TA 
exp (- Teh 
Yo 




















x1— 


oe. D) ei (23) 


oC 


in which (23) is proposed to approximate (22). The tightness of this approximation will be verified later in the 
next section. Herein, g denotes the mean interference channel power gain at the CUE and is calculated as 

















Dit+tRe l 
g Ê Ell] = / ai 24) 
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5. NUMERICAL RESULTS 

In this section, Monte Carlo simulations are constrcuted to validate our analysis in the previous sec- 
tions. The default settings for the results are Re = 15 (m), Ra = 5 (m), Rin = 0.5 (bps/Hz), % = 10 (dB) 
and ņ = 2. In addition, the average transmit SNR is set to 7 = 30 (dB) in the Figure 6. 
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Figure 5. Outage probability at the CUE versus the average transmit signal-to-noise of the BS (7) and the 
distance D; (m) 


Figure 5 shows the outage performance comparison of conventional user as changing D+. It can be 
seen that as varying D, from 5 (m) to 12.5 (m), a small performance gap emerges. Increasing transmit SNR 
at the BS contributes to improving outage performance. It is confirmed that analytical and simulation results 
are matched very tight at the whole range of SNR. Hence, it can be believed that the proposed approximation 
acceptably tight. Figure 6 further evaluates the mean interference channel power gain of the considered system 
to validate our exact analysis for the link distance distribution. 
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Figure 6. Mean interference channel power gain at the CUE versus the distance D; (m) 





CONCLUSION 
In this paper, we have studied the outage performance problem of a cellular user in the context of 


the D2D underlying cellular network. The outage probability of users in D2D communication is computed. 
Specifically, the distance from the cellular user to the BS and from the D2D transmitter to the user are both 
statistically characterized using stochastic geometry tools. We also derived the exact and approximated expres- 


sions 


of outage probability at the cellular node. Numerical results are provided to verify our derivations. The 


paper is, however, limited at downlink single-cell networks, and will be extended to the multi-cell sencario in 
our future works. 
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